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ANIMAL CAGE BEHAVIOR SYSTEM 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims priority to and benefit of USSN 60/399,043, filed on 

July 25, 2002, which is incorporated herein by reference in its entirety for all purposes. 

STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER FEDERALLY 
SPONSORED RESEARCH AND DEVELOPMENT 
[0002] This invention was made with Government support by Grant Nos: MH61624 

and MH01949 awarded by the National Institutes of Health. The Government of the United 
States of America has certain rights in this invention. 

FIELD OF THE INVENTION 
[0003] This invention pertains to the field of ethology and neurobiology. In 

particular this invention provides devices and methods to evaluate genetic influences on 
animal behavior. 

BACKGROUND OF THE INVENTION 

[0004] Substantial limitations currently exist in our ability to apply recent 

biotechnological advances to analyze neural substrates of complex mammalian behavior. In 
contrast to the rapid pace of innovation seen in the fields of mammalian genomics, 
medicinal chemistry and information technology, less progress has been made in the 
development of behavioral assessment techniques for mice or other mammals. Such 
procedures are vital for exploring the impact of genes, drugs and environment on brain 
functions relevant to common neuropsychiatric conditions such as schizophrenia, 
depression, and anxiety. Standard approaches involving repeated removal of mice from 
their home cages for a battery of behavioral tests are problematic because: 1) they are time- 
consuming and labor-intensive, 2) the order of test administration can skew the resulting 
data, 3) removal of mice from the home cage produces stress that confounds interpretation 
of behavioral data, and 4) data are frequently misinterpreted due to a failure to consider 
behavioral domains that are not the main focus of study (eg: impact of anxiety on tests of 
learning). 
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[0005] For example, the 2-bottle ethanol preference test (often used in behavioral 

phenotyping in alcohol research) is widely used in genetic screens to provide an indication 
of the rewarding properties of ethanol. The test is simple to administer, but difficult to 
interpret; increased preference could either reflect enhanced sensitivity to the rewarding 
properties of ethanol, or alternatively, a reduction of sensitivity requiring animals to 
consume more ethanol to achieve similar hedonic effects. Moreover, preference for 
ethanol-containing solutions can be influenced by factors such as taste, novelty and caloric 
content, and it is often unclear whether animals achieve blood alcohol levels sufficient to 
produce pharmacological effects. Studies dependent on isolated behavioral outcome 
measures such as preference are therefore subject to overinterpretation if they are not placed 
in the context of the multiple neuropsychological processes that influence behavioral 
performance. 

[0006] The problem of reliance on isolated behavioral readouts also pertains to 

apparently simple measures, as indicated by recent findings that the impact of genetic 
mutations on ethanol-induced motor impairments are highly dependent on the particular 
behavioral assays employed. In humans, alcohol abuse, dependence and withdrawal occur 
in a complex behavioral context and produce global perturbations of behavioral 
organization, including alterations in sleep and vigilance, motor activity levels and 
coordination, nutrition, anxiety and mood. Behavioral phenotyping approaches are needed 
that address the complex behavioral contexts in which alcohol consumption occurs, as well 
as those behaviors that are in turn influenced by alcohol. 

SUMMARY OF THE INVENTION 
[0007] This invention pertains to a novel behavioral monitoring system useful for 

the analysis of complex behaviors in a number of animal species. The monitoring systems 
of this invention allow continuous monitoring of feeding, drinking and movement of 
animals with high temporal and spatial resolution. The high temporal and spatial resolution 
the systems of this invention permits the detection, quantification, and categorization of 
various behavior patterns, and/or components of such patterns and facilitates automated 
analysis and generation of behavioral fingerprints. 
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DEFINITIONS 

[0008] The term "lickometer" refers to a device that provides a measure of fluid 

consumption by an animal. 

[0009] The term "feeder" refers to a device that provides food to an animal in 

captivity. In preferred embodiments, the feeder provides a measure of the amount of food 
consumed by the animal. 

[0010] The term "database" refers to a means for recording and retrieving 

information. In preferred embodiments the database also provides means for sorting and/or 
searching the stored information. The database can comprise any convenient media 
including, but not limited to, paper systems, card systems, mechanical systems, electronic 
systems, optical systems, magnetic systems or combinations thereof. Preferred databases 
include electronic (e.g. computer-based) databases. Computer systems for use in storage 
and manipulation of databases are well known to those of skill in the art and include, but are 
not limited to "personal computer systems", mainframe systems, distributed nodes on an 
inter- or intra-net, data or databases stored in specialized hardware (e.g. in microchips), and 
the like. 

[0011] A "behavioral fingerprint" refers to a summary of behavioral components 

that is characteristic of an animal under certain defined conditions and that can be used to 
distinguish that animal from another species and/or strain and/or mutant. 

[0012] An "event file" is a file containing behavioral data where each entry is 

characterized by an event. This is in contrast to a "continuous recording file". The event 
file can comprise any medium capable of recording information, e.g., paper, optical media, 
magnetic media, internet storage, etc. In preferred embodiments, the event file comprises a 
computer readable medium. 

[0013] The term "test agent" refers to an agent that is to be screened in one or more 

of the assays described herein. The agent can be virtually any chemical compound. It can 
exist as a single isolated compound or can be a member of a chemical (e.g. combinatorial) 
library. In a certain embodiments, the test agent will be a small organic molecule. 

[0014] The term "small organic molecule" refers to a molecule of a size comparable 

to those organic molecules generally used in pharmaceuticals. The term excludes biological 
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macromolecules {e.g., proteins, nucleic acids, etc.). Preferred small organic molecules 
range in size up to about 5000 Da, more preferably up to 2000 Da, and most preferably up 
to about 1000 Da. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 [0015] Figure 1 schematically illustrates a behavior analysis system 40 comprising 

an enclosure 38, an activity platform 20, a lickometer 32 comprising a water bottle 26 and a 
detector 24, a feeder 18 comprising a ramp 12, a detector 10 {e.g. a photobeam), food 14, 
and an overflow trough 16. The system is illustrated with a test animal 22 {e.g. a mouse), 
and an optional video and/or image acquisition system 28. In this illustration, the enclosure 

10 38 is topped with a cover 30. In certain embodiments, feeding is measured by a photobeam 
and photobeam detector with 1 msec resolution when animals break the photobeam in order 
to reach food in the feeder. Animals accessed the feeder via a ramp 12. Licking is 
measured at the lickometer with, e.g., a 1 msec resolution by a change in capacitance in the 
licking spout when the animal licked to obtain water. In certain embodiments, the feeder 

15 and lickometer are on the same side of the cage. Movement can be measured on the activity 
platform by a load beam detector capable of detecting movement of greater than 1 cm with 
20 msec resolution. 

[0016] Figure 2 shows a photograph of a behavior analysis system 40. 

[0017] Figure 3 shows a schematic illustration of the feeder 18 showing the entrance 

20 ramp 12, the detector 10 (e.g. photobeam), the feed bin 42, and the overflow trough 16. 

[0018] Figures 4A and 4B show the correlation of chow intake with beam break 

time for 5-HT 2C mutant (knockout mice (Fig. 4A), and for wild type mice (Fig, 4B). 

[0019] Figure 5 shows a schematic illustration of the lickometer 32 showing the 

fluid reservoir 26, the detector 24 (e.g. conoductance or capacitance meter), and the barrier 
25 plate 52, with slot 54. 

[0020] Figures 6A and 6B show the correlation of fluid intake with number of licks 

5-HT 2C wild type mutant (knockout mice (Fig. 6A), and for 5-HT 2C mutant (knockout 
mice (Fig, 6B). 
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[0021] Figure 7 illustrates feeding, drinking, and movement of an animal in a 6 

minute interval. 

[0022] Figure 8 illustrates a representation of the spatial organization of mouse 

behavior monitoring device according to this invention. The graph shows the duration of 
5 pauses between movement events for one dark cycle displayed based on location within the 
home cage. 

[0023] Figure 9 illustrates three consecutive N-N excursions. 

[0024] Figure 10 illustrates an analysis of activity type. 

[0025] Figure 1 1 illustrates the structure of a feeding bout or meal. 

10 [0026] Figure 12 illustrates an exponential fits to frequency of pauses between 

movements for different types of movements. 

[0027] Figure 13 shows behavioral fingerprints (ethograms) for various mouse 

species. 

[0028] Figure 14 shows behavioral fingerprints (ethograms) for various mutant 

15 (knockout) mice. 

[0029] Figure 15 shows behavioral fingerprints (ethograms) for various mutant 

(knockout) mice. 

[0030] Figure 16 shows a photograph of a behaioral monitoring system comprising 

multiple feeders and multiple lickometers. 

20 [0031] Figure 17 shows correlations of intake with device measures. Left plot: 

Photobeam break time was highly correlated with chow intake measured in grams for wild 
type (R 2 =.9606) and 2C mice (R 2 =.9631). Right plot: Lick contact number was highly 
correlated with water intake measured in grams for wild type (R 2 =.9602) and 2C mice 
(R 2 =.9845). Measurements were done for 3, 6, 9, 12, and 24 hour periods and across time 

25 of day. 

[0032] Figures 18A, 18B, and 18C show total consumption and movement in the 

light and dark cycles. Chow and water intake and amount of movement in the light and 
dark cycles were analyzed using repeated measures ANOVA to allow for comparisons of 
the effects of genotype, cycle, and cycle by genotype interactions. Figure 18A: Total chow 
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consumed per cycle. 2C, Ob, and Ay mice had significant effects of cycle (p<.001; p<.001; 
p<.001) as well as significant effects of genotype (p<.01; p<.01; p<.001) on the amount of 
chow consumed when compared to WT mice. 2C mice had a significant cycle by genotype 
interaction (p<.05). Figure 18B: Total water consumed per cycle. 2C, Ob, and Ay mice 
5 had significant effects of cycle on the amount of water consumed (p<.001; p<.001; p<.001). 
2C and Ob mice had significant cycle by genotype interactions (p<.01;p<.01). Figure 18C: 
Total movement per cycle. 2C, Ob, and Ay mice had significant effects of cycle (p<.001; 
p<.001; p<.001), significant effects of genotype (p<.01; p<.001; p<.01) and had significant 
cycle by genotype interactions (p<.03; p<.001; p<.001) on total movement per cycle. In all 
10 figures, values represent mean ± SEM, with wild type indicated by grey symbols and 
mutant indicated by black symbols. 

[0033] Figure 19 twenty-four hour feeding pattern for wild type, 5HT2CR, Ob, and 

Ay mutant mice. The light cycle occurs between hours 0 to 12 and the dark cycle occurs 
between hours 12 and 24, and a dotted line represents the transition from light to dark cycle. 
15 A red line represents each feeding event. Bout and cluster, identified using their respective 
criteria, are labeled in green and purple respectively. 

[0034] Figure 20 shows frequency histograms of photo beam break gaps on a scale 

of log minutes with fitted log normal distributions. Histograms are of single mice with 9 to 
10 days of data and are separated by dark and light cycles. Individual log normal 
20 distributions are shown in blue and the best fit is shown in red. The intersections of log 

normal distributions are marked with green dotted lines. The cluster criteria are assigned at 
the intersection of the last and second to last distributions, and the bout criteria are assigned 
at the intersection of the second to last and third to last distributions. The intersections of 
distributions are marked with green dotted lines. 

25 [0035] Figure 21 shows the probability of being in a position during a gap. During 

within-bout intervals (WBI), animals had the highest probability of being at the feeder. 
During inter-bout intervals (IB I), animals had equal probabilities of being at the feeder, at 
the licking spout, and engaged in activity around the cage. During inter-cluster intervals 
(ICI), animals had the highest probability of being at the nest. The x- and y-axes are 

30 measured in cm. The z-axis is the probability of being at a position, with the scale shown 
on the colorbar. The square represents the feeder and the circle represents the licking spout. 
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[0036] Figures 22A, 22B, and 22C shows cluster properties separated by cycle and 

genotype. Cluster properties were analyzed using T-tests to compare each mutant group 
with their respective wild type controls for each cycle. Figure 22A: Mean cluster size per 
cycle, measured by amount of food consumed in grams. In the dark cycle, Ob mice had 
5 significantly larger cluster sizes (p<.01) than their WT controls. In the dark cycle, Ob and 
Ay mice had significantly larger cluster sizes (p<.001;p<.001) than their respective WT 
controls. Figure 22B: Mean number of clusters per cycle. In the dark cycle, 2 C and Ob 
mice had significantly fewer clusters (p<.05; p<.02) than their respective WT controls. In 
the light cycle, 2C mice had significantly more clusters (p<.01) and Ob mice had 
10 significantly fewer clusters (p<.01) than their respective WT controls. Figure 22C: Mean 
intercluster interval duration per cycle. In the dark cycle, Ob mice had a significantly 
longer ICI duration (p<.001) than their WT controls. In the light cycle, 2C mice had 
significantly shorter ICI durations (p<.001) and Ob mice had significantly longer ICI 
durations (p<.001) than their respective WT controls. 

15 [0037] Figure 23 shows cluster size measured in grams of food consumed in 6 

minute bins across circadian time. The light cycle occurs between hours 0 to 12 and the 
dark cycle occurs between hours 12 and 24. Circadian rhythms in ingestive behavior can be 
seen by the variations in cluster sizes between the dark and light cycles and ultradian 
rhythms can be seen by the variations in cluster sizes in the dark cycle. Each red square 

20 represents a cluster of feeding. 

[0038] Figures 24A and 24B shows periodicities in feeding behavior of wild type 

versus Ob mutant mice. Figure 24A; The amounts of food eaten in 6 minute intervals for 
ten days of data were averaged and smoothed, shown in red. Figure 24B: Periodograms, 
showing the frequencies of certain periodicities of feeding across time, were calculated from 
25 the averaged and smoothed data. The dotted lines show the significances values. The 

significant peaks were used to generate a waveform, shown in blue applied to the data in 
(Figure 24A. WT mice have significant peaks at 24, 12, 8, and 5 hours. Ob mice have 
significant peaks at 24 and 3 hours. 

[0039] Figure 25A, 25B, and 25C show clustering based on periodograms and 

30 averaged waveforms generated from significant periods. Figure 25A: Kullback-Leibler 

measures of divergence were calculated between periodograms for individual mice and used 
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to cluster based on Ward hierarchical clustering technique. The clustering of periodograms 
for individual mice shows three qualitative groups, one with mostly WT and 2C mice, one 
with mostly WT and Ay mice, and one with only Ob mice. Figure 25B: The averaged 
periodograms for each group are shown. Figure 25C: Waveforms generated from each 
5 averaged periodogram are shown. 

[0040] Figure 26 illustrates a representational digital device 700 according to the 

present invention. 

DETAILED DESCRIPTION 

L Behavioral monitoring system. 

10 [0041] This invention pertains to a novel behavioral monitoring system useful for 

the analysis of complex behaviors in a number of animal species. The monitoring systems 
of this invention allow continuous monitoring of feeding, drinking and movement of 
animals with high temporal and spatial resolution. The high temporal and spatial resolution 
the systems of this invention permits the detection, quantification, and categorization of 

15 various behavior patterns, and/or components of such patterns. 

[0042] The components of the system (e.g. the lickometer, the feeder, etc.) are 

optimized with a number of features that that minimize or significantly reduce the incidence 
of erroneous readings. The devices therefore provide an output that is a good measure of 
actual animal behavior. 

20 [0043] Thus, for example, in certain embodiments, the detector in the feeder 

comprises a photobeam. It was a surprising discovery that small animals (e.g. rats, mice, 
etc.) tend to pile up food during a feeding activity and this piled up food can block the 
photobeam, or other detector and produce a false positive reading. Such false positives 
were reduced/eliminated by engineering the feeder with a trough 16 around the feed bin 42 

25 (see, e.g., Figure 3). As illustrated by Figures 4A and 4B, detector signal (e.g. photobeam 
breaks) provided a good measure of actual food consumption. 

[0044] The lickometer 32 is also modified to reduce false positives. As illustrated in 

Figure 5, in certain embodiments, the lickometer is shielded from the animal by a plate 52 
comprising an opening or slot 54. The plate prevents the animal from simply touching the 
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lickometer, e.g. with its paws. The signals produced by lickometer are therefore simply a 
function of the animal licking the spout rather than other manipulations of the device. As 
illustrated by Figures 6A and 6B, detector signal (e.g. licks) provided a good measure of 
actual fluid consumption. 

5 [0045] In various embodiments, animal position and/or activity is monitored using a 

motion or position detection system. Suitable motion or position detection systems include, 
but are not limited to photobeam arrays, pressure transducers, force plates, video and/or 
photographic image systems, and the like. In certain embodiments, activity monitoring is 
accomplished using a two dimensional force plate to provide location in X and Y 
10 coordinates. 

[0046] The monitoring systems of this invention operate at a significantly higher 

spatial resolution and temporal resolution than typical monitoring systems. It was a 
surprising discovery that such high resolution data acquisition/data reporting provides a 
result different in kind than that observed with typical systems. Specifically, the high 
15 temporal and spatial resolution permits the detection and recording of behavioral 

components that could not previously be detected and/or quantified. Moreover, the high 
temporal and spatial resolution permits curve fitting and quantititative statistical analysis of 
resulting data sets permitting automated and/or semiautomated detection, quantification, 
and/or characterization of behavioral components. 

20 [0047] In certain embodiments the monitoring systems of this invention report (log) 

data with a temporal resolution of greater than 20 seconds, preferably greater than 10 or 1 
second, more preferably greater than 100 ms, or 50 ms, and most preferably greater than 25 
ms, 10 ms, or 1 ms. 

[0048] Spatial resolution is typically better than 10 cm, preferably better than about 

25 5 cm, more preferably better than about 3 cm or about 2 cm, and most preferably better than 
about 1 cm, 0.5 cm, or 0.1 cm. 

[0049] The systems of this invention allow the investigation of multiple behavioral 

components while leaving the test animals relatively undisturbed. This affords a number of 
advantages over previous testing schemes. In particular, the systems of this invention allow 
30 the investigation of behaviors in context and the simultaneous investigation of multiple 

behavioral domains. The systems of this invention avoid the effects of test order on animal 
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behavior observed with test batteries. The animals are not subject to the stress of handling 
and removal from their home cage. The freezing of animals (e.g. observed in zero maze 
tests) is not observed in these systems. 

II, Data collection and creation of an ethogram. 
5 [0050] As indicated above, the systems of this invention provide behavioral data and 

high spatial and temporal resolution. To prevent "data swamping", in certain embodiments, 
data is not recoded continuously, but rather at the occurrence of events to produce an "event 
file". An event file is record delineating a collection of events. 

[0051] Thus, for example, a position/locomotion event file will record an event 

10 every time the animal moves a distance greater than a certain predetermined distance (e.g. 1 
cm). Each time the animal moves a distance greater than the predetermined distance, the 
event file can record, e.g. the animal location and the time or the duration from the previous 
event. 

[0052] Similarly, for drinking or eating behavior, the event file can record the time 

15 of duration between signals from feeder or lickometer. 

[0053] The data thus collected can be used to provide a detailed portrayal of the 

animals complex behavior during any particular time period (see, e.g. y Figure 7). Thus, for 
example, Figure 7 shows illustrates feeding, drinking, and movement of an animal in a 6 
minute interval.. Figure 8 illustrates a representation of the spatial organization of mouse 
20 behavior monitoring device according to this invention. The graph shows the duration of 
pauses between movement events for one dark cycle displayed based on location within the 
home cage. Figure 9 illustrates three N-N excursions, and Figure 10 illustrates a specific 
activity analysis. 

[0054] The animal(s) behavior can be broken down into bouts and clusters where 

25 bouts are repeated occurrences of the same act (e.g., photobeam breaks or food 

consumption) that appear to cluster together in time (see, e.g., Figure 1 1). Clusters are 
collections of bouts that can define a particular behavior (e.g. feeding, grooming, etc.). The 
intervals separating the repeating acts give the bout structure and the intervals separating the 
bouts give the cluster structure. Thus, for example, feeding bouts or "meals" have been 
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described as having two types of intervals: Within-meal intervals and between-meal 
intervals. 

[0055] The quantitative data can be analyzed using various numeric and/or 

statistical methods (see, e.g., Figure 12) to provided detailed information regarding 
5 behavioral patterns. Moreover, this analysis can readily be automated. 

[0056] The collection of bouts and clusters can be analyzed to provide a behavioral 

fingerprint of the animal at a particular time and/or under particular conditions. 

[0057] It was a surprising discovery that a characteristic behavioral finger print can 

be determined for a particular animal strain or species, or particular knockout animals, or 
10 for animals subject to agents (e.g. pharmaceuticals) that alter or potentially alter behavioral 
patterns. 

[0058] Thus, for example, Figure 13 shows behavioral fingerprints (ethograms) for 

various mouse species, while Figures 14 and 15 show behavioral fingerprints for various 
mousse mutants. 

15 [0059] Data can be collected from the position detector(s), and/or lickometer, and/or 

feeder using standard instrumentation well known to those of skill in the art. As indicated, 
herein, in certain preferred embodiments, the feeder uses breakage of a photobeam to detect 
feeding, while the lickometer uses changes in capacitance. Circuits for detecting 
photobeam breakage and/or capacitance changes are well known to those of skill in the art 

20 and are commercially available. Thus, for example, the ENV-254 "head entry" package 

produced by MED Associates, Inc. (St. Albans, Vermont ) is designed to detect photobeam 
breaks by intrusion of a mouse head into a receptacle. MED Associates also provides a 
variety of photobeam and/or capacitance circuits for use in lickometers. The various units 
work with controllers provided by the same company that readily interface with 

25 commercially available data acquisition systems. These units are illustrative and not 

intended to be limiting. Those of skill in the art can readily design custom circuitry for use 
in the devices of this invention. In addition a number of other detectors are commercially 
available. 
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A) Analytical approaches for behavioral pattern analysis. 

1. Classification of behavioral elements and their aggregation into 
clusters 

[0060] One ethologically-based approach is to develop criteria for defining 

5 behavioral bouts, where bouts are defined as repeated occurrences of the same behavioral 
act, which are themselves temporally aggregated into clusters (Machlis (1977) Behaviour 1- 
70). In one approach, a mixture of probability density distributions are fit to a large set of 
inter-event intervals using maximum likelihood estimation {see, e.g., Example 2). The 
feeding data that has been analyzed at this time indicate that the data for most mice are best 

10 fit by three to six inter-event intervals. Typically, the first one to three intervals are brief 
(less than 30 seconds), and seem likely to reflect within-bout intervals where the animal 
continues to feed and does not leave the feeder. In addition, videotape observation of mice 
feeding at the feeder suggests that the different brief inter-event intervals seen during a bout 
of feeding may reflect differences in how the mouse handles the food (using its paws to 

15 remove the food and eat versus eating directly with its mouth while its head is in the food 
drawer). The next to last distribution typically has a mean of around one minute and seems 
likely to represent the inter-bout interval. The presence of the last distribution indicates that 
the feeding bouts themselves are clustered together in time. This approach is readily 
automated to apply bout and cluster criteria to the behavioral record (e.g. as illustrated in 

20 Example 2). 

[0061] It was a surprising discovery that this technique can be used to classify the 

behavior of the mouse into bouts of feeding, drinking, and movement and in addition, a 
higher level of organization — clusters of bouts, can also be detected. The preliminary work 
was done by visually inspecting an inter-event interval frequency histogram for each mouse 

25 and making initial estimates of the parameters of the fit that were then used to initiate an 

iterative fitting routine (modified Newton-Raphson) (Press et al. (1992) Numerical Recipes 
in C, Cambridge University Press, New York). This was necessary because fitting 
algorithms like Newton-Raphson are sensitive to the initial estimate provided, because 
while they work well to identify local maxima, they often fail to find the global maxima. 

30 Also, maximum likelihood estimation provides no model selection criteria, i.e. there is no 
robust analytical way to choose the optimal number of distributions to use in the mixture. 
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In order to solve this problem and automate this procedure, we use a technique developed 
by Maneesh Sahani for the sorting of electrophysiological data (Sahani (1999) Latent 
variable models for neural data analysis. Thesis, California Institute of Technology). This 
method has a number of benefits relative to other algorithms—it implements a global 
5 maximization algorithm similar in spirit to methods such as a deterministic annealing 
expectation maximization, and split and merge expectation maximization. Unlike these 
methods, however, it allows simultaneous automatic selection of the best number of 
components in the mixture (clusters). 

[0062] The strategy may then be used to identify bouts of not only feeding but also 

10 of drinking and movement in an analogous automated high throughput manner. 



2. Characterization of circadian influences on cluster properties 
[0063] Inspection of home cage behavioral patterns reveals that behavioral clusters 

are strongly modulated by circadian factors. For example, cluster durations are briefer 

15 during the light cycle than the dark cycle for C57BL/6 mice (see, e.g., Example 2). 

Moreover, even within the dark cycle, cluster sizes of C57BL/6 mice vary in a systematic 
way~they are largest early and late in the dark cycle. Circadian influences on behavioral 
patterns vary greatly among strains, many of which do not display the large light- 
onset/light-offset clusters observed in the C57BL/6 strain. For example, A/J mice do not 

20 appear to have large dark-offset clusters, and 129Sl/SvImJ animals do not form large 
clusters at any time of day. Single gene energy balance mutations also influence the 
periodicity of feeding, as shown by loss of the C57BL/6 pattern and replacement by clusters 
displaying less diurnal variability in size, as well as an apparent 3-hr periodicity. 

[0064] We thus observe that behavioral clusters vary with time, exhibit both 

25 circadian (24 hr periodicity) and ultradian (less than 24 hr periodicity) rhythms, and are 
highly sensitive to genetic influences. These important features of the behavioral patterns 
can be quantified using the methods described herein.. In one illustrative approach, data is 
binned and converted to the frequency domain using the Lomb-Scargle algorithm. This 
approach, developed for the analysis of astrophysical data, has several advantages (Ruf 
30 (1999) Biol. Rhythm Res. 30: 178-201). These include the capacity of the algorithm to deal 
with missing data points. Also, the Lomb-Scargle algorithm provides a robust method for 
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determining the statistical significance of the peaks found in the periodograms 
(periodograms provides visual depictions of the extent to which behavioral events are 
organized into periods of various durations). These benefits make this algorithm a good 
choice power spectrum estimator for this data set. 

5 3. Characterization of within-cluster behavioral patterns 

[0065] In addition to objectively defining behavioral clusters and the diurnal 

modulation of their properties, the methods described herein permit examination of patterns 
of behavioral units within clusters. One approach is to hierarchically assess the transition 
probabilities between successively longer sequences of behavior. For example, movement, 

10 feeding and drinking can be designated as three distinct behavioral states. Because, in 

certain embodiments, the feeder and lickometer are separated, at different sites in the cage, a 
feeding animal has a higher probability of transitioning to a movement state than to drinking 
(since the mouse must move to reach the lickometer). Progressing from this somewhat 
trivial example, the analysis can be extended consider the relative probabilities of 

15 progressively longer sequences of behavior. For example, it is possible to determine 
whether the feeding^movement^feeding sequence is more or less probable than the 
feeding— >movement^drinking sequence. 

[0066] In another embodiment, an alternative approach for identifying behavioral 

patterns that utilizes concepts from information theory, a computational approach enabling 
20 quantitative assessment of the extent to which events are organized. Elements of patterned 
behaviors are more predictive of future actions than those of random or disorganized 
behaviors. It is possible to therefore calculate the predictive information for recorded data 
sequences. This allows quantitative assessment of the extent to which experimental 
manipulations alter the patterning of spontaneous behavior. 

25 [0067] The characterization of within-cluster behavior can also be further clarified 

by a more detailed analysis of movement. First, an attempt will be made to classify 
locomotor behavior with regard to the speed and turning angle of the movements. The 
focus here would be to identify movement states such as walking, running, or lingering. In 
one approach, pause durations between movements will be fit with log normal distributions 

30 and assessments performed of the extent to which groups of movements generated from the 
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different distributions have distinct properties such as differences in spatial spread or 
turning angle. Another strategy involves application of a predictive information approach to 
a series of movements in a search for abrupt changes in information content reflecting 
changes in the movement type. 

5 [0068] Spatial information can also be used to examine the behavior of the mice 

during intervals between bouts of ingestive behavior and/or sleep. In certain "ethograms", 
bouts of feeding are separated by interbout intervals (IB Is) of different lengths. The - 
locomotor paths taken by this animal during the IBIs (3 round-trips from the feeder) are 
distinct. For example, in certain ethograms, the first reflects a round-trip between the feeder 
10 and lickometer. The second incorporates a visit to the next. The third reflects a different 

behavior that tends to occur near the end of bouts-animals sweep bedding along the outside 
of the niche, gathering material that will be used to seal the opening to the niche. Once 
derived, these paths can be treated as images and clustered using standard pattern 
recognition algorithms. 

15 4. Classification algorithms for group comparisons 

[0069] Once the behavioral elements and their organization have been characterized 

for individual mice, comparisons can be made between treatment groups. The large volume 
and complexity of behavioral data and the inability to control all of the determinants of 
behavioral variability create challenges in distinguishing uncontrolled individual differences 

20 from group effects arising from planned experimental manipulations. Several statistical 

modeling techniques can be applied to this problem. One approach utilizes hidden Markov 
models (HMMs) (see, e.g., Rabiner (1989) Proc IEEE 77: 257-285). Hidden Markov 
models are useful for modeling time-dependant data. The heart of a hidden Markov model 
is statistical in nature — the underlying system generating the data is modeled as a statistical 

25 process with a number of "hidden" states that make transitions from one to another in such a 
way that the probability that the system is in one state at a certain time is only dependant on 
what state the system was in at the previous time step. These hidden states then generate 
the observable sequence of behavioral units, such as feeding, drinking and movement bouts. 
To generate the model, various optimization techniques are used to find the HMM that best 

30 reflects the data. To generate the model, one can use standard techniques such as the Baum- 
Welch algorithm to find the HMM that best reflects the data. These outputs provide a much 
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more compact representation of the data and can be used to quantify the extent to which 
behavioral patterns differ. The ability to make such distance ("relatedness") measurements 
permits the use of standard clustering techniques (hierarchical clustering, self-organizing 
maps, etc) to identify groups of animals with similar patterns of behavior. Post-hoc 
5 analyses of the particular features of behavioral patterns that distinguish groups can provide 
insights into the neurobehavioral processes that are altered by experimental manipulations. 

[0070] While the foregoing analytical approaches are described in detail with 

respect to mice, it will be appreciated that, using the teaching provided herein, they can 
readily be applied to other mammals {e.g. rats, cats, dogs, primates, etc.). 

10 B) Embodiment in a programmed digital apparatus 

[0071] The invention may be embodied in a fixed media or transmissible program 

component containing logic instructions and/or data that when loaded into an appropriately 
configured computing device cause that device to perform one or more of the analytical 
operations described above on an ethological dataset {e.g. classify behavior into bouts, 

15 identifiy circadian patterns to behavioral bouts, classivy within cluster behaviors, compare 
groups, etc.) according to the methods of this invention. 

[0072] Figure 26 shows digital device 700 that may be understood as a logical 

apparatus that can read instructions from media 717 and/or network port 719. Apparatus 
700 can thereafter use those instructions to direct analysis of behavioral data, create, sort, 
20 search, and read bevaioral database, and the like. In certain embodiments, the digital device 
can be directly connected to one or more cage behavioral systems according to this 
invention and, optionally function in realtime. In certain embodiments, the digitial device 
can simply access, analyze, and/or manipulate previously collected data. 

[0073] One type of logical apparatus that may embody the invention is a computer 

25 system as illustrated in 700, containing CPU 707, optional input devices 709 and 711, disk 
drives 715 and optional monitor 705. Fixed media 717 can be used to program such a 
system and could can represent disk-type optical and/or magnetic media, and/or a memory 
or the like. Communication port 719 can also be used to program such a system and can 
represent any type of communication connection {e.g. a connection to a data acquisition 
30 system). 
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[0074] The invention also may be embodied within the circuitry of an application 

specific integrated circuit (ASIC) or a programmable logic device (PLD). In such a case, 
the invention may be embodied in a computer understandable descriptor language that can 
be used to create an ASIC or PLD that operates as herein described. 

5 C. Embodiment in a web site. 

[0075] The methods of this invention can be implemented in a localized or 

distributed computing environment. In a distributed environment, the methods can be 
implemented on a single computer comprising multiple processors or on a multiplicity of 
computers. The computers can be linked, e.g. through a common bus, but more preferably 
10 the computer(s) are nodes on a network. The network can be a generalized or a dedicated 
local or wide-area network and, in certain preferred embodiments, the computers may be 
components of an intra-net or an internet. 

[0076] In certain internet embodiments, a client system typically executes a Web 

browser and is coupled to a server computer executing a Web server. The Web browser is 

15 typically a program such as Microsoft's Internet Explorer, or NetScape or Opera. The Web 
server is typically, but not necessarily, a program such as IBM's HTTP Daemon or other 
WWW daemon. The client computer can be bi-directionally coupled with the server 
computer over a line or via a wireless system. In turn, the server computer can be bi- 
directionally coupled with a website (server hosting the website) providing access to 

20 software implementing the methods of this invention. 

[0077] A user of a client connected to the Intranet or Internet can cause the client to 

request resources that are part of the web site(s) hosting the application(s) providing an 
implementation of the methods of this invention. Server program(s) then process the 
request to return the specified resources (assuming they are currently available). A standard 

25 naming convention has been adopted, known as a Uniform Resource Locator ("URL"). 
This convention encompasses several types of location names, presently including 
subclasses such as Hypertext Transport Protocol ("http"), File Transport Protocol ("ftp"), 
gopher, and Wide Area Information Service ("WAIS"). When a resource is downloaded, it 
may include the URLs of additional resources. Thus, the user of the client can easily learn 

30 of the existence of new resources that he or she had not specifically requested. 
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[0078] The software implementing the method(s) of this invention can run locally 

on a server hosting the website in a true client-server architecture. Thus, the client 
computer posts requests to the host server which runs the requested process(es) locally and 
then downloads the results back to the client. Alternatively, the methods of this invention 
can be implemented in a "multi-tier" format wherein a component of the method(s) are 
performed locally by the client. This can be implemented by software downloaded from the 
server on request by the client (e.g. a Java application) or it can be implemented by software 
"permanently" installed on the client. 

[0079] In one embodiment the application(s) implementing the methods of this 

invention are divided into frames. In this paradigm, it is helpful to view an application not 
so much as a collection of features or functionality but, instead, as a collection of discrete 
frames or views. A typical application, for instance, generally includes a set of menu items, 
each of with invokes a particular frame-that is, a form which manifest certain functionality 
of the application. With this perspective, an application is viewed not as a monolithic body 
of code but as a collection of applets, or bundles of functionality. In this manner from 
within a browser, a user would select a Web page link which would, in turn, invoke a 
particular frame of the application {i.e., subapplication). Thus, for example, one or more 
frames may provide functionality for inputing and/or accessing ethograms for particular 
animals or strains, , while another frame provides tools for identifying bouts, clusters, 
circadian patterns, and the like. 

[0080] In addition to expressing an application as a collection of frames, an 

application can also be expressed as a location on the Intranet and/or Internet; a URL 
(Universal Resource Locator) address pointing the application. Each URL preferably 
includes two characteristics: content data for the URL (i.e., whatever data is stored on the 
server) together with a data type or MIME (Multipurpose Internet Mail Extension) type. 
The data type allows a Web browser to determine how it should interpret data received from 
a server (e.g., such as interpreting a .gif file as a bitmap image). In effect, this serves as a 
description of what to do with the data once it is received at the browser. If a stream of 
binary data is received as type HTML, the browser renders it as an HTML page. If instead it 
is received type bitmap, on the other hand, the browser renders it as a bitmap image, and so 
forth. 
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[0081] In Microsoft Windows, different techniques exist for allowing a host 

application to register an interest in a data object {i.e., data of a particular type). One 
technique is for the application to register with Windows an interest in a particular file 
extension for an {e.g., .doc-"Microsoft Word Document"); this is the most common 
5 technique employed by Window applications. Another approach, employed in Microsoft 
Object Linking and Embedded (OLE), is the use of a class Globally Unique Identifier or 
GUK)~a 16-byte identifier for indicating a particular server application to invoke (for 
hosting the document having the GUID). The class ID is registered on a particular machine 
as being connected to a particular DLL (Dynamic Link Library) or application server. 

10 [0082] In one embodiment of particular interest, a technique for associating a host 

application with a document is through a use of MIME types. MIME provides a 
standardized technique for packaging a document object. It includes a MIME header for 
indicating which application is appropriate for hosting the document, all contained in a 
format suitable for transmission across the Internet. 

15 [0083] In one preferred embodiment, the methods of the present invention are 

implemented, in part, with the use of a MIME type specific to the use of the methods of this 
invention. The MIME type contains information necessary to create a document (e.g., 
Microsoft ActiveX Document) locally but, in addition, also includes information necessary 
to find and download the program code for rendering the view of the document, if 

20 necessary. If the program code is already present locally, it need only be downloaded for 
purpose of updating the local copy. This defines a new document type which includes 
information supporting downloadable program code for rendering a view of the document. 

[0084] The MIME type may be associated with a file extension of . APP. A file with 

the .APP extension is an OLE Document, implemented by an OLE DocObject. Because the 

25 .APP file is a file, it can be placed on a server and linked to using an HTML HREF. The 
.APP file preferably contains the following pieces of data: (1) the CLSID of an ActiveX 
object, which is an OLE Document Viewer implemented as one or more forms appropriate 
to the use of the methods of this invention; (2) the URL of the codebase where the object's 
code can be found, and (3) (optionally) a requested version number. Once the .APP 

30 DocObject handler code is installed and registers the APP MIME type, it can be used to 
download an .APP file into the user's Web browser. 
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[0085] On the server side, since the . APP file is really a file, the Web server simply 

receives the request and returns the file to the client. When the APP file is downloaded, the 
.APP DocObject handler asks the operating system to download the codebase for the object 
specified in the .APP file. This system functionality is available in Windows through the 
5 CoGetClassObjectFromURL function. After the ActiveX object's codebase is downloaded, 
the .APP DocObject handler asks the browser to create a view on itself, for instance, by 
calling the ActivateMe method on the Explorer document site. The Internet Explorer then 
calls the DocObject back to instantiate a view, which it does by creating an instance of the 
ActiveX view object from the code that was downloaded. Once created, the ActiveX view 
10 object gets in-place activated in the Internet Explorer, which creates the appropriate form 
and all its child controls. 

[0086] Once the form is created, it can establish connections back to any remote 

server objects it needs to perform its functions. At this point, the user can interact with the 
form, which will appear embedded in the Internet Explorer frame. When the user changes to 
15 a different page, the browser assumes responsibility for eventually closing and destroying 
the form (and relinquishing any outstanding connections to the remote servers). 

[0087] In one preferred embodiment, from an end-user's desktop, the entry point to 

the system is the corporate home or the home page of another particular web-site. The page 
can, optionally, include, in a conventional manner, a number of links. In response to the 
20 user clicking on a particular link to an application page {e.g. a page providing the 

functionality of the methods of this invention), the web browser connects to the application 
page (file) residing on the server. 

[0088] In one embodiment, where the user requests access to the methods of this 

invention, the user is directed to a particular page type, e.g., an application (appdoc) page 

25 for in-place execution of an application (implementing one or more elements of the methods 
of this invention) in the Web browser. Since each application page is located using an 
URL, other pages can have hyperlinks to it. Multiple application pages can be grouped 
together by making a catalog page that contains hyperlinks to the application pages. When 
the user selects a hyperlink that points to to an application page, the Web browser 

30 downloads the application code and executes the page inside the browser. 
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[0089] Upon the browser downloading the application page, the browser (based on 

the defined MIME type) invokes a local handler, a handler for documents of a type, ore 
particularly, the application page preferably includes a Globally Unique Identifier (GUID) 
and a codebase URL for identifying a remote (downloadable) application to invoke for 
5 hosting the document. Given the document object and the GUID which arrive with the 
application page, the local handler looks to the client machine to see if the hosting 
application already resides locally (e.g., by examining Windows 95/NT registry). At this 
point the local handler can choose to invoke a local copy (if any) or download the latest 
version of the host application. 

10 [0090] Different models of downloading code are commonly available. When code 

is downloaded, a "code base" specification (file) is initially requested from the server. The 
code base itself can range from a simple DLL file to a Cabinet file (Microsoft .cab file) 
containing multiple compressed files. Still further, an information (e.g., Microsoft .inf) file 
can be employed for instructing the client system how to install the downloaded application. 

15 These mechanisms afford great flexibility in choosing which component of an application 
gets downloaded and when. 

[0091] For preferred embodiments, the machinery employed for actually 

downloading program code itself relies on standard Microsoft ActiveX API (Application 
Programming Interface)-calls. Although the ActiveX API does not provide native support 

20 for Web-delivered applications, its API can be invoked for locating the correct version of 
the program code, copying it to the local machine, verifying its integrity, and registering it 
with the clients operating system. Once the code has been downloaded, the handler can 
proceed to invoke the now-present application host for rendering the document object (in a 
manner similar to invoking the hosting application through the registry if it were already 

25 installed). 

[0092] Once the hosting application (OLE server) is loaded at the client, the client 

system can employ the OLE document view architecture to render the application correctly 
within the browser, including using conventional OLE methodology for adding the 
application's menu to that of the browser and for correctly re-sizing the application upon a 
30 re-size of the browser (as oppose to requiring the application to execute within a single 
Active X control rectangle—the limitation previously noted). Once the application is 
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executing at the client, it can execute remote logic such as using RPC (Remote Procedure 
Call) methodology. In this manner logic which is preferably implemented as remote 
procedure(s) can still be used. 

[0093] In certain preferred embodiments, the methods of this invention are 

5 implemented as one or more frames providing the following functionality. Function(s) to 
organize, search, save, and retrieve raw behavioral data or reduced/processed behavioral 
data {e.g. data produced by the devices of this invention), functions to identify and/or 
classify bouts, functions to identify/classify clusters of bouts, functions to identify/classify 
circadian patterns, functions to classify/identify within bout behaviors, functions to compare 
10 and contrast ethograms, functions to graphically represent ethograms, and the like. 

[0094] In addition, the functions can also, optionally, provides access to private 

and/or public databases accessible through a local network and/or the intranet whereby one 
or more ethograms contained in the databases can be input into the methods of this 
invention. 

15 [0095] Methods of implementing Intranet and/or Intranet embodiments of 

computational and/or data access processes are well known to those of skill in the art and 
are documentede in great detail {see, e.g., Cluer et al. (1992) A General Framework for the 
Optimization of Object-Oriented Queries, Proc SIGMOD International Conference on 
Management of Data, San Diego, California, Jun. 2-5, 1992, SIGMOD Record, vol. 21, 

20 Issue 2, Jun., 1992; Stonebraker, M., Editor; ACM Press, pp. 383-392; ISO-ANSI, 

Working Draft, "Information Technology-Database Language SQL", Jim Melton, Editor, 
International Organization for Standardization and American National Standards Institute, 
Jul. 1992; Microsoft Corporation, "ODBC 2.0 Programmer's Reference and SDK Guide. 
The Microsoft Open Database Standard for Microsoft Windows.TM. and Windows 

25 NT.TM., Microsoft Open Database Connectivity.TM. Software Development Kit", 1992, 
1993, 1994 Microsoft Press, pp. 3-30 and 41-56; ISO Working Draft, "Database Language 
SQL-Part 2:Foundation (SQIVFoundation)", CD9075-2:199.chi.SQL, Sep. 11, 1997, and the 
like). 

[0096] Those skilled in the art will recognize many modifications can be made to 

30 this configuration without departing from the scope of the present invention. For example, 
in a two-tier configuration, the server system executing the functions of the WWW gateway 
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may also execute the functions of the Web server. For example, any one of the above 
described embodiments could be modified to accept requests from users/user terminals that 
are in a format other than a URL. Yet another modification would involve the adaptation to 
a multi -manager environment. 

5 HI. Uses of the bahavioral monitoring system and ethograms. 

[0097] The equipment and methods described herein can be used in a wide variety 

of contexts. For example, the creation of ethograms for various knockout animals provides a 
direct link between genetic variations and behavioral repertoire. Moreover, the combination 
of such analyses of knockout animals in conjunction with the administration of various 
10 compounds, provides a powerful mechanism for dissecting the biological basis of 

overconsumption, addiction, and the like and for identifying the genetic component(s) 
associated with particular behavioral patterns. 

[0098] In addition, the methods of this invention can readily be used to investigate 

possible psychoactive properties of essentially any test agent. Such an evaluation can 
15 become a critical component of any regulatory approval process for pharmaceuticals, food 
additives, compounds released into the environment, and the like. 

[0099] Such methods can, for example, involve contacting (administering to) a 

animal subject with one or more test agents of interest. One or more behavioral finger 
prints (ethograms) under defined conditions can then be determined for the animal and 
20 compared to appropriate controls (e.g. ethgrams determined for a similar negative control 
animal (e.g. an animal administered less or no test agent, or various positive controls). 

[0100] The controls can be determined at the time the test experiment is done, 

and/or cam be ethograms done previously and possibly compiled into a reference ethogram 
database. 

25 [0101] The following illustrates one use of the equipment and methods of this 

invention: 

[0102] Quantitative approaches to behavioral pattern assessment can be used to 

examine the impact of ethanol on the organization of home cage behavior in the mouse. 
One series of studies examines the behavioral consequences of a standard 2-bottle voluntary 
30 ethanol intake situation. Ethanol consumption patterns and their effects on home cage 
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behavioral organization are correlated with blood ethanol concentrations. To assess levels 
of intoxication, standard tests of ethanol-related behaviors cam be performed at times of 
peak blood ethanol concentrations. A subsequent analogous series of studies can be 
performed in animals rendered tolerant to ethanol by prior exposure to ethanol vapor. 
5 Following stabilization of behavioral patterns, effects of ethanol withdrawal on home cage 
behavioral organization can be determined and correlated with standard behavioral indices 
of withdrawal. These tests can be applied to mouse strains known to exhibit high and low 
levels of ethanol preference. 

[0103] C57BL/6J and DBA/2 J strains can be used to determine how the effects of 

10 ethanol on behavioral patterns compare in high and low ethanol-pref erring strains that have 
been used extensively in genetic studies of alcohol preference (46). Although C57BL/6J 
mice display high ethanol preference, it is unclear whether they consume ethanol for its 
pharmacological properties or for its taste or caloric content (47). To determine whether 
mice are drinking alcohol for its pharmacological effects, it is necessary not only to measure 

15 the patterns of intake, but also to know the resultant concentrations of alcohol that are 
achieved in the blood. We can thus apply a metabolic model for predicting the blood 
ethanol concentrations (BEC) of a mouse based on measurement of its rates of intake, 
distribution, and elimination (4). C57BL/6J and DBA/2J mice will be given 10% ethanol 
by intra-gastric gavage, and blood samples collected to estimate the rates of absorption, 

20 distribution, and elimination for use in calibrating the metabolic model. The capacity to 
continuously predict the BEC of mice using continuous monitoring of intake can then be 
tested. Mice are placed in the behavioral monitors in the presence of a water bottle and a 
10% ethanol bottle. Once intake is stable, mice are removed from the monitors at several 
times of day BEC determinations, and the degree of correlation of predicted levels with 

25 ethanol concentrations measured from blood samples will be tested. Corrections to the 

model will be made at this point if necessary, and the experiment repeated until an adequate 
capacity to predict BEC from the rate of ingestion of 10% ethanol using continuous 
monitoring has been achieved. 

[0104] Because genetic studies of ethanol intake frequently involve 2-bottle testing 

30 using 10% ethanol vs water (48), we can examine the impact of ethanol on home cage 
behavior using this paradigm. C57BL/6J and DB A/2J mice will be exposed to this 
situation, along with control groups that receive water in both bottles. Bout and cluster 
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analysis of ethanol intake can be examined along with analysis of the influence of ethanol 
on home cage behavioral patterns. We can use our data on the ingestion of ethanol to 
continuously predict the BEC of each mouse and compare these predicted levels across 
strains and experimental day. These BEC level predictions can be plotted against activity, 
5 drinking, and feeding for initial visual examination of differences in the patterns of ethanol 
intake, and its impact on other behaviors across the strains. We can correlate patterns of 
movement, derived from the behavioral analysis, with BEC in an attempt to detect signs of 
intoxication. At the end of the data collection period, subsets of mice from each strain are 
removed from their cages at times predicted to correspond to maximal levels of BEC. 

10 Animals can be tested on routine measures of intoxication including rotorod testing of 
coordination and elevated zero-maze testing of anxiety levels. After these tests, blood 
samples will be taken and compared to the values predicted by monitoring of consumption. 
The performance of test mice and control mice removed from cages at the same times will 
then be compared to assess whether or not the test mice from each strain voluntarily 

15 consume sufficient quantities of alcohol to exhibit significant signs of intoxication. 

[0105] It is not certain that voluntary consumption of the 10% ethanol 

solution will result in sufficient intake to produce intoxication and dependence. To examine 
the impact of ethanol dependence and withdrawal on behavioral organization, mice can be 
rendered ethanol-dependent using chronic exposure to ethanol vapor (CEEV) (49, 50). 

20 C57BL/6J and DBA/2J mice can be placed in cages housed within chambers equipped to 
deliver ethanol vapor (test group) or water vapor (control) for four weeks. BEC can be 
periodically monitored to assess the adequacy of vapor exposure (49-52). Monitoring will 
can be carried out over 2 weeks, and we expect to observe significant increases in voluntary 
ethanol intake in the chronic ethanol vapor exposed group. Alterations in the patterns of 

25 ingestion responsible for this increase in intake can be determined, and the impact of this 
increased intake on patterns of home cage behavior can be analyzed. Testing for 
coordination impairment and reduction of anxiety during predicted peak BEC levels can be 
performed in one group of mice. Ethanol can be removed from the lickometers of another 
group of animals to determine behavioral patterns associated with ethanol withdrawal. We 

30 would anticipate reductions of sleep episodes, anorexia, adipsia and fragmentation of 
activity patterns (53, 54). 
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[0106] Another example of the use of the devices of this invention is the analysis of 

the impact of serotonin receptor mutations on ethanol intake, dependence and ethanol influences on 
behavioral organization. Although brain serotonin systems strongly influence ethanol-related 
behaviors, the contributions of particular 5-HT receptor subtypes to these actions of serotonin are 
5 unclear. We have generated 3 lines of null mutant mice, each lacking a 5-HT receptor subtype that 
may significantly influence ethanol responses. The devices of this invention can be used to evaluate 
the impact of these mutations on behavioral patterns associated with ethanol intake and dependence. 
The lines to be tested are: 1) 5-HT2CR mutants, that exhibit enhanced mesolimbic dopaminergic 
activity and elevated cocaine self-administration, 2) 5-HT1AR mutants, that display a global 

10 elevations of serotonergic tone and enhanced anxiety-related behavior and 3) 5-HT6R mutants, that 
display reduced sensitivity to acute behavioral effects of ethanol and enhanced ethanol self- 
administration. Although changes in ethanol pharmacokinetics are not anticipated in these mutants, 
this will be confirmed experimentally. Separate cohorts of mutant and wild type mice will then be 
tested in the standard 2-bottle voluntary ethanol intake situation with and without prolonged 

15 exposure to ethanol vapor in the manner described in Aim 2. 

EXAMPLES 

[0107] The following examples are offered to illustrate, but not to limit the claimed 

invention. 

Example 1 

20 [0108] In certain embodiments, test animals (e.g. mice, rats, rabbits, guineae pigs, 

etc.) in cages equipped with a feeding monitor, a lickometer, and an activity detection 
platform (see, e.g. A Figures 1 and 2). To determine food and fluid intake, feeding and 
licking coefficients (FC, LC) are calculated by dividing total consumption by the total 
photobeam break (PBB) time and lick contact duration (LCD). Multiplication of the FC by 

25 individual PBB times then yields the instantaneous food consumption while multiplication 
of the LC by LCD yields the instantaneous water consumption. We have validated this 
approach, observing tight correlations between measured intakes and PBB and LCD times. 

[0109] This system has been used to characterize home cage behavioral patterns in 

C57BL/6J, hyperphagic ob/ob (26) and 5-HT2CR null mutant mice (27, 28). Inspection of 
30 a record from a C57BL/6J mouse reveals that feeding is episodic, closely associated with 

both water consumption and movement, and expressed in a diurnal pattern. Examination of 
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a record from an ob/ob mouse reveals less frequent, but more intense episodes of ingestive 
behavior, consistent with impairment in satiation processes. In addition, a marked decrease 
in movement is apparent; a factor likely to contribute to the enhanced adiposity of ob/ob 
mice. Interestingly, the pattern of hyperphagia in 5-HT2CR mutant mice differs from that 
5 of ob/ob mice, characterized by increased frequencies of ingestive bouts, raising the 

possibility that the regulation of hunger (processes leading to feeding initiation) is perturbed 
in these animals. 

[0110] We have demonstrated that home cage behavior for mice is organized and 

sensitive to mutations that influence neural function. Detailed analysis of behavioral 
10 organization provides a highly sensitive readout of the influences of genetic, 
pharmacological and environmental manipulations on CNS function. 

[0111] We have developed quantitative ethological approaches to understand the 

organization of behavior in the mouse. The development of these basic techniques provides 
a means to determine patterns of alcohol ingestion, or ingesion of other compounds, in a 
15 broad behavioral context, and to examine how this ingestion affects behavioral organization, 
and by extension, brain function. 

[0112] Episodes of ingestive behavior are characterized by complex clusters 

of eating and drinking bouts interspersed with locomotion and other behaviors (29-31). To 
perform pattern analysis of consummatory behaviors, objective criteria are developed for 

20 defining these bouts and clusters. Toward this end, we have used maximum likelihood 
estimation to fit a log-normal model to the frequency distribution of intervals between 
feeding events (31-33). An example of the distribution of the inter-event-intervals for 
feeding for an individual mouse reveal that a best fit for the feeding events utilizes five 
gaussians. Our movement data indicate that mice do not leave the feeder during the brief 

25 intervals comprising the first three distributions of the feeding event data. Intervals between 
24 seconds and 18 minutes fall within the fourth distribution, and constitute intervals 
between bouts of feeding. Frequently, these intervals correspond to visits to the licking 
spout during which drinking episodes ensue. The longest pauses of greater than 18 minutes 
define intervals between clusters of feeding and drinking bouts (Fig 3C) and frequently 

30 corresponds to visits to the nest. A similar approach may be used to uncover structure in 
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and locomotor behavior to provide information relating to the size, duration, frequency and 
diurnal distribution of bouts and clusters of behavior. 

[0113] Another important organizing feature of home cage behavior is the nest. 

Wild mice establish nests to which they repeatedly return after excursions into their home 
5 range (34, 35). Our preliminary data indicate that singly housed laboratory mice use the 
nest as a home base in a similar manner. The longest pauses between movements cluster in 
one area of each cage that corresponds to the nest (Fig 4A). We examine individual 
excursions from this location by defining nest-to-nest loops which describe the location of 
the mouse and its behavior during a round-trip from the nest and back. We observe that 
10 these trips frequently correspond with visits to the feeder and the water bottle where feeding 
and drinking bouts ensue. 

[0114] This inevntion provides quantitative behavioral analysis techniques to define 

patterns of movement that provide sensitive indicator(s) of ethanol intoxication or the 
influence of other psychoactive chemicals. Ethological studies have described multiple 

15 activity states, including full stop, lingering, walking, turning and running (35, 36). We 

have used a method similar to that described by Drai et al to define states of movement that 
we predict will correspond to running, walking or turning, moving in place, and full arrest. 
The example of a nest-to-nest excursion reveals a pattern we commonly observe, in which 
velocity of movement away from the nest is less than the velocity exhibited during the 

20 return trip. 

5-HTR Mutant Lines. 

[0115] Consistent with a contribution of 5-HT2CRs to the tonic serotonergic 

inhibition of mesolimbic dopamine system function, 5-HT2CR mutant mice display 
elevated extracellular DA levels in the nucleus accumbens, as determined using a no-net 

25 flux in vivo microdialysis procedure. Several lines of evidence suggest an association of 

mesoaccumbens DA projection hyperactivity with disinhibited exploration of environmental 
novelty (37, 38). Accordingly, 5-HT2CR mutants exhibit disinhibited exploration of a 
moderately aversive mirrored chamber and elevated locomotor activity in a novel open 
field. Enhanced responses to novelty have been associated with increased sensitivity to the 

30 locomotor and reinforcing effects of psychostimulants (39). Accordingly, 5-HT2CR 

mutants exhibit increased locomotor responses to cocaine (Fig 5D) and will lever press at 
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twice the rate of wild type mice for IV cocaine infusions in an operant progressive ratio 
task. In view of these studies, we predict that 5-HT2CR activation opposes the reinforcing 
properties of ethanol. 

[0116] Consistent with an inhibitory role of 5-HT1A somatodendritic autoreceptors 

5 on brain serotonergic activity, 5-HT1 AR mutants exhibit global elevations of extracellular 
serotonin levels (40). Consistent with Cloninger's model relating hyperserotonergic states to 
anxiety, these animals also display a profound anxiogenic phenotype (41). 

[0117] Unlike the 5-HT2CR and 5-HT1 AR mutant lines, no overt baseline 

behavioral abnormalities were detectable in 5-HT6 null mutants. However, these animals 

10 display markedly abnormal responses to acute ethanol administration, characterized by 
reduced ataxic (rotorod assay) and sedative (loss of righting reflex) effects, and enhanced 
locomotor stimulation . These mutants overconsume ethanol solutions in chronic intake 
studies, achieving elevated blood ethanol levels. Neural mechanisms underlying this 
phenotype are unclear; the relatively high levels of 5-HT6R expression in the nucleus 

15 accumbens (23, 24) provides a focus for further study. 

Example 2 



Introduction 

[0118] The production of behavior by an animal does not occur randomly in time 

20 nor space, but is instead organized in the context of the environment to improve the fitness 
and thus the survival of the animal (Machlis (1997) Behaviour 63: 1-70). As a result of this 
organization, discernable and coherent patterns of behavior should be detectable, however 
the continual production of behavior by an animal in the natural or laboratory environment 
produces a vast amount of information. The challenge therefore is not only to develop the 
25 technology to obtain the behavioral data, but also to develop analytical techniques that 
reduce the information contained in this behavioral data to a succinct and informative 
summary. Because of the role of evolutionary selection in behavioral adaptation (Tinbergen 
(1963) 20: 89-94), there should be a substantial genetic contribution to the central nervous 
system structures that produce behavioral patterns. A detailed description and 
30 quantification of such behavioral patterns should therefore elucidate the genetic and 
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neurophysiological mechanisms that regulate the generation of behavior. As a step in 
developing automated high-throughput techniques for the collection and analysis of 
behavioral patterns, a home cage behavioral monitoring system for mice that captures with 
high temporal and spatial resolution the feeding, drinking, and locomotion of an animal has 
5 been developed. 

[0119] Mice are good experimental subjects for these studies because of the 

sequencing of their genome, the availability of inbred strains and spontaneous mutations, as 
well as the increasing ability to manipulate their genetics. The rationale for examining 
feeding, drinking, and locomotion results from the prediction that activities essential for 

10 survival can be expected to strongly influence the organization of behavior. The 

maintenance of energy balance is essential for survival, and depends on the regulation of 
energy intake through alterations in feeding and drinking behavior and on energy 
expenditure, which is affected behaviorally by levels of physical activity. The role of 
feeding, drinking, and physical activity levels in the maintenance of energy balance thus 

15 suggests that these behaviors should be highly organized in time and space. In addition, the 
availability of mutations producing obesity provides an excellent opportunity to evaluate the 
impact of alterations in energy balance on behavioral patterns. 

[0120] Three mutations {obese, agouti, and 5HT2c receptor) that disrupt the 

regulation of energy balance and lead to obesity via alterations in leptin, melanocortin, and 

20 serotonergic systems (Pellymounter et al (1995) Science 269: 540-543; Fan et al (1997) 
Nature 385: 119-120; Tecott et al (1995) Nature 374: 542-546) have been investigated 
using the home cage behavioral monitoring system. Of these, the obese mutation results in 
marked obesity due to the loss of function of the hormone leptin that is expressed primarily 
in adipose tissue and acts at receptors in the ventrobasal hypothalamus as well as other sites 

25 (Boston et al (1997) Science 278: 1641-1644; Schwartz et al (1996) Journal of Clinical 
Investigation 96: 1101-1 106). In the absence of leptin, obese mice increase feeding and 
decrease activity (Pellymounter et al (1995) Science 269: 540-543). In addition, they are 
hyperglycemic, hyperinsulinemic, hypothermic, and have an overall decrease in metabolic 
rate (Friedman and Halaas (1998) Nature 395: 763-770), which contributes to the profound 

30 obesity seen in obese mice. The agouti mutation also produces marked obesity but results 
from ectopic expression of the agouti protein that acts as an antagonist at the melanocortin 
type 3 and 4 receptors (Miller et al (1993) Genes and Development 7: 454-467; Oilman et 
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al (1997) Science 278: 135-138). These receptors are highly expressed in the 
paraventricular and dorsomedial hypothalamic nucleus and in the lateral hypothalamic area 
(Kim et al (2000) Diabetes 49: 177-182). Obesity in double mutant agouti/obese mice has 
been shown to be independent and additive suggesting that the mechanism leading to 
5 obesity in agouti mice is distinct from that leading to obesity in obese mice (Boston et al. 
(1997) Science 278: 1641-1644). Finally, the engineered null mutation of the serotonin 2c 
(5HT2C) receptor results in mild late onset obesity, accompanied by development of type II 
diabetes. The 5HT2C receptor is widely expressed in the central nervous system but has 
high levels of expression in the arcuate nucleus of the hypothalamus and the nucleus 

10 accumbens. While all three mutations lead to increased food intake and obesity, the 
mechanisms leading to obesity appear to be distinct, suggesting that the alterations in 
ingestive behavior and physical activity levels produced by these mutations may also be 
distinct. Analysis of patterns of feeding, drinking, and physical activity combined with 
comparison of patterns between different obese mutants thus provides a means to assess the 

15 utility and adequacy of the description of behavioral patterns. 

[0121] In order to identify patterns of behavior exhibited by an animal, it is first 

useful to reduce the observation of the continuous flow of behavior into simple behavioral 
units. A common unit used to measure behavior in ethology is the bout, which is defined 
as the repeated occurrence of the same behavioral act clustered together in time (Machlis 

20 (1997) Behaviour 63: 1-70). The underlying notion behind this definition is that the animal 
should be engaged in a single behavior during a bout (e.g. feeding) and that termination of 
this behavior should be separated by the animal engaging in other distinct behaviors (e.g. 
drinking, movement, or inactivity). In addition, individual behavioral acts within a bout, 
such as repeatedly grasping food with the paw and bringing it to the mouth, are separated by 

25 short within-bout intervals (WBIs) that are distinct in duration from the longer intervals that 
separate bouts. It is the existence of both short WBIs and longer inter-bout intervals (IBIs) 
that produces the clustering in time that is characteristic of a bout, thus if bouts of behavior 
are present, it should be possible to demonstrate that at least two types of inter-event 
intervals are present separating the repetition of a behavior. This can be done by fitting a 

30 mixture of probability density distributions to a large set of inter-event intervals using 
maximum likelihood estimation to assess the presence of distinct types of inter-event 
intervals (Tolkamp et al (1998) Journal of Theoretical Biology 194: 235-250). If the 
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frequency of inter-event intervals is better fit by a mixture of two distributions this indicates 
that two distinct durations of intervals are present and supports the proposal that the 
behavior is organized into bouts. In addition, the fit to the data can also be used to 
determine a bout criteria- the maximum duration of time that an animal can pause between 
5 the repeated behavioral acts before the bout is considered to have ended. This criterion can 
then be used to identify the bouts of behavior and thus classify the continuous stream of 
behavior into behavioral units. The bout is a useful unit of classification because 
presumably the initiation and termination of bouts is regulated by internal and external cues 
to which the animal responds. 

10 [0122] Two types of mixture distributions (exponential and log normal) have 

typically been fit to the frequency distribution of inter-event intervals when attempting to 
identify the presence of a bouts in the production of a behavior. In the past, exponential 
distributions have frequently been fit to the data, however the use of exponentials requires 
the assumption that the probability of initiating a feeding event is independent of the time 

15 since the last feeding event. Furthermore, the use of exponentials implies that the 

probability of initiating feeding is highest directly after a bout, and decreases exponentially 
thereafter. This is not consistent with the idea of satiety, which would suggest that the 
animal has a low probability of initiating feeding immediately following a bout of feeding. 
Tolkamp and Kyriazakis (1999) Animal Behaviour 57: 807-817, have proposed fitting the 

20 inter-event intervals with log normal distributions, which means that the probability of 

initiating feeding following a bout is low and increases with time until reaching the mode of 
the log normal distribution. This model is therefore a more meaningful model as it is 
consistent with the concept of satiety. 

[0123] In addition to identifying the occurrence of bouts in feeding behavior, these 

25 types of analysis also revealed that bouts of feeding may themselves be clustered together in 
time indicating that a higher level of organization is present (i.e. clusters of bouts). The 
presence of clusters of bouts is indicated when the frequency distribution of inter-event 
intervals is best fit by a mixture of three distributions. This corresponds to the presence of 
short duration WBIs, intermediate duration EBIs, and longer duration inter-cluster intervals 
30 (ICIs). These clusters are often called meals, and variation in the properties of the clusters, 
such as their duration and the quantity of food consumed during the cluster, is related to 
biologically phenomena such as satiety/hunger and satiation. 
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[0124] As a first step in examining the patterns of feeding behavior in the obese 

mutants examined in these experiments, log normal distributions are fit to the intervals 
between feeding acts to assess for the presence of bouts and clusters. The properties of 
these bouts and clusters are then compared across genotypes to examine how these 
5 mutations alter the ingestive behavior patterns of these obese mutant mice. The role of 
circadian and ultradian rhythms in determining the patterns of ingestion is also examined. 

Methods 

Animals 

[0125] In experiments involving 5HT2C receptor mutant mice, hemizygous mutant 

10 males bearing a null mutation of the X-linked htrlc gene (congenic on a C57BL6J 

background) were used with littermate wild types (WT) as controls. This line has been 
maintained by mating females heterozygous for the htrlc gene with C57BL6J males 
(Jackson Laboratories, Bar Harbor, ME) and has been genotyped as previously reported 
(Tecott et al (1995) Nature 374: 542-546). In experiments involving obese mice, 

15 homozygous mutant mice (B6.V -Lep ob 9 Jackson Laboratories, Bar Harbor, ME) were used 
with WT mice from the Jackson Laboratory C57BL/6J stock as controls. In experiments 
using agouti mice, heterozygous mutant mice (B6.Cg-A y , Jackson Laboratories, Bar Harbor, 
ME) were used with WT from the Jackson Laboratory B6.Cg-A y stock as controls. All 
mice were two-three months of age when the studies began. Animals were always 

20 maintained on a 12-hour light/dark cycle (lights off at 1900 h) and were given free access to 
water and standard chow diet (PicoLab Mouse Diet 20; Purina Mills, Richmond, IN). The 
average body weights of Ob (57.0g±3.3) and Ay (40.2g±.0.9) mice during the course of the 
experiment differed significantly (p<.001;p<.001) from those of their wild types 
(28.1g±1.9). The average body weights of 5HT2CR mutant mice were not significantly 

25 different than those of their wild type controls. Animals were used and cared for according 
to the University of California, San Francisco laboratory animal research guidelines. 

Data Collection 

[0126] During each experiment, mice were singly housed in 20 cm by 50 cm cages 

with Cell Sorb Plus Litter Bedding (Fangman Specialties Inc., Cinncinati, OH) and a 
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bedding square. Cages were mounted on a load beam platform (DiLog Instruments, 
Tallahassee, FL) that detected movement of the animal's center of gravity greater than 1 cm 
with a 20 ms temporal resolution. Each cage had a feeder and a licking bottle (DiLog 
Instruments, Tallahassee, FL) attached to one end (see, e.g., Figure 1). The feeder allowed 
5 access to a powered version of the standard chow (PicoLab Mouse Diet 20; Purina Mills, 
Richmond, IN) through a hole in the wire mesh ramp leading into the feeder. Feeding was 
detected when an animal broke an infrared photobeam mounted just below the hole in the 
feeder ramp. Drinking was detected as a change in capacitance when an animal licked at the 
lick spout of the water bottle. The temporal resolution of the photobeam break and lick 

10 contact measurements was 1 ms. The validity of using photobeam break time and lick 

contact number as a measure of food and water intake was assessed by measuring food and 
water intake after varying periods of time (3, 6, 9, 12, and 24 hours). Total intake and the 
device measures of feeding and licking were highly correlated for wild type (R 2 = 9606; 
R 2 =.9602) and 5HT2CR (R 2 =.9631; R 2 =.9845) mutant mice across all times of day (Figure 

15 17). Photobeam beak time and lick contact number thus provide a reasonable measure of 
food and water intake. 

[0127] Mice were run in groups of 16 with 8 mutants (either 5HT2CR null, obese, 

or agouti) and 8 controls (2 5HT2CR rounds, 1 Ob round, and 1 Ay round). An additional 
round of 9 Ob mice was run with heterozygote controls which were excluded from analysis 

20 due to their partial phenotype compared with wild type controls. Each round lasted 14 days 
consisting of 4 days for habituation and 10 days for data collection. Data was collected 
continuously from the devices except during daily maintenance during which food and 
water were removed and replaced with fresh food and water. Food and water were 
weighted before and after removal to measure total intake for the duration of data collection. 

25 At the end of one round using 5HT2CR mutants and wild types, correlation tests were done 
to assess the validity of the devices for measuring intake as described above. 



Data Analysis 

Quality Control 

[0128] The accuracy of the device measurements was assessed at the end of each 

30 round to determine if errors in measurement were present. The most problematic error in 
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data collection stemmed from the powdered food either piling up or sticking to the 
photobeam LED or IR detector resulting in prolonged photobeam breaks unrelated to 
feeding. In order to exclude these erroneous photobeam breaks, the photobeam break data 
was binned into 30 second bins and all bins that had 30 seconds of continuous beam break 
5 time and at least one centimeter of locomotion were combined and the frequency of 
movement was fit with two exponentials. A criterion was then set as the intersection 
between the two exponentials. It was assumed that movement less than the criterion resulted 
from movement in the feeder during the photobeam break and that movement greater than 
the criterion resulted from movement around the cage during an erroneous photobeam 

10 break. Therefore if an animal moved more than the criterion (9-1 1cm depending on the 

round) during a 30 second bin of continuous beam break time, the entire day of data for that 
animal was excluded. This technique appears effective since all data containing a 
simultaneous photobeam break and lick contact were detected and removed from the data 
set. An additional problem occurred more rarely and resulted when an animal slept in the 

15 feeder. This was first determined by visually assessing the feeding records for abnormal 

amounts of photobeam break time and then observing the animal in the day. The entire day 
of data was excluded if the animal was observed to sleep in the feeder. Additionally, one 
day of data for an entire round of Ob mice was lost due to a temperature malfunction that 
caused a prolonged temperature fluctuation to 85°F. Finally, an animal was completely 

20 dropped from analysis if, as a result of the exclusions, it had four or fewer days of data 
available for analysis. 

Estimation of Intake 
[0129] Intake for any given time duration less than that measured directly was 

calculated using feeding and licking coefficients. The feeding coefficient (mg/s) was 

25 calculated by dividing the total intake for the duration of data collection (usually about 22 
hours) by the total photobeam break time for the duration. The licking coefficient (mg/lick) 
was calculated by dividing the total intake for the duration of data collection by the total 
lick contact number for the duration. Multiplying the feeding coefficient by the total 
photobeam break time (sees.) during the duration of interest (e.g. dark or light cycle) then 

30 yields an estimate of the quantity of food consumed during that duration. A similar 
calculation yields an estimate of the quantity of water consumed. 
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Bout and Cluster Analysis 
[0130] All photobeam break gaps (intervals between photobeam breaks) were 

compiled for each animal for the last 10 days of data collection. Due to obvious differences 
in the patterns of the gaps between the light and dark cycle, the gaps for the dark and light 
5 cycle for each animal were examined separately. The distribution of the dark and light 
cycle log-transformed gaps for each animal were then individually fit with one to nine log 
normal distributions by minimizing the deviance (deviance = -2*X(0)) of the log-likelihood 
(A.(0), below) for a mixture of log normal distributions using the modified Newton-Raphson 
method in GENSTAT: 

10 Me) = Zm-nZm-mC ln(pj* exp(-(t-m) 2 /(2*a j 2 ))/(27ra j 2 ) 1/2 

where pj is the proportion of intervals in distribution j, |Jj is the mean of distribution j, Gj is 
the standard deviation of distribution j, t is the log(interval duration), with interval duration 
in milliseconds, N is the number of intervals, and M is the number of distributions. 
[0131] The number of distributions fit to the data was increased until the maximum 

15 likelihood ratio, LR = 2*(^(0i)-X(e 2 )), fell below 15. The last fit with a ratio above 15 was 
deemed significant and accepted as the final fit. 

[0132] The distribution with the longest mean interval duration was the distribution 

of inter-cluster intervals (ICI). The distribution with the second longest mean interval 
duration was the distribution of inter-bout intervals (IB I). The remaining distributions, with 

20 the shortest mean interval durations, were collectively termed the distributions of within- 
bout intervals (WBI). The criteria for assignment of interval type were determined by the 
intersections of distributions. The cluster criterion was calculated as the intersection 
between the IBI distribution and the ICI distribution. Intervals with duration longer than the 
cluster criterion were labeled ICIs. The bout criterion was calculated as the intersection of 

25 the IBI with the adjacent intersecting WBI, which is not always the next shortest duration 
due to variation in the standard deviations. Intervals with a duration longer than the bout 
criterion but shorter than the cluster criterion were IB Is and those with duration shorter than 
the bout criterion were labeled WBIs. Separate bout and cluster criterion were calculated 
for the dark and light cycle. 

30 [0133] Some exceptions to this method of classifying the gaps were made to provide 

for consistent comparison across animals. For some animals (n = 2, 2C, light cycle), in the 
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region where the IBI distribution was typically found (-0.5-2 minutes), there were two 
distributions. This suggests, that for these animals, there were two types of IBI or that the 
IBI changed over the 10 days of data collection. In this case, the bout criterion was 
recalculated at the intersection of the shorter IBI and the adjacent intersecting WBI. 
5 Animals were also said to have no cluster criteria if the highest duration intersection of 
curves was below 5 minutes (n = 5, 2C, dark cycle). These animals either lacked clusters 
during the dark cycle or there were too few ICIs during the dark cycle and the ICI was 
therefore not detected. However, even if there were no clusters during the dark cycle, all of 
these animals did exhibit an ICI distribution during the light cycle. The dark cycle could 

10 thus be considered as a single large cluster separated from clusters occurring during the 
light cycle by the light cycle ICIs. A criterion of 720 minutes (12 hours of the dark cycle) 
was therefore compared with a criterion yielded by averaging the criteria of animals of the 
same genotype. The criterion of 720 minutes would be consistent with no distinct clusters 
appearing during the dark cycle, while the application of an average criterion from other 

15 mice of the same genotype is consistent with too few ICIs being present for detection. 

Inspection of the data revealed that application of the averaged criterion identified clusters 
that seemed consistent with those estimated by visual inspection. It thus seemed more 
likely that for these animals cluster did occur during the dark cycle but that there were too 
few ICIs to be detected. Therefore, an averaged criterion was used for animals without 

20 cluster criteria. 

[0134] After the interval types were assigned, the photobeam events themselves 

were assigned as belonging to bouts or clusters, depending on which types of intervals 
surrounded them. This assignment was done for all photobeam breaks using Visual Basic 
for applications in Microsoft Excel. For analysis, bouts and clusters were then separated 
25 into active and inactive phases, because clusters sometimes crossed the boundary of the 

dark-light cycle. Clusters that began in the light cycle but ended in the dark cycle or began 
in the dark cycle but ended in the light cycle were assigned to the active phase (largely dark 
cycle). Clusters that occurred solely within the light cycle were assigned to the inactive 
phase (light cycle). 
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Time Series Analysis 

[0135] Total intake during 6 minute bins was calculated for 10 days of data. The 

Lomb-Scargle algorithm was used to generate a periodogram for each animal and 
significant peaks from each periodogram were used to generate a waveform. The 
5 divergence between periodograms was measured using the Kullback-Leibler measure of 
distance and the periodograms were clustered based on these distances using the Ward 
hierarchical clustering technique. 

Statistical analysis 

[0136] Total food and water intake and total movement in the light and dark cycles 

10 were compared using repeated measures ANOVA. Cluster properties, including cluster size 
(grams), number of clusters, and inter-cluster interval, were compared using repeated 
measures ANOVA. Statistical tests were run separately for 2C, Ob, and Ay mice and 
mutants were compared to the control WT mice for each group. 

Results 

15 Alterations in Diurnal Patterns of Ingestive Behavior and Activity 

[0137] Because mice are known to be nocturnal, differences in the temporal 

organization of behavior between the dark and light cycles were expected. Therefore, total 
food and water consumption as well as movement in the dark and light cycles were first 
examined. Not surprisingly for all groups of animals, there was a significant effect of cycle 

20 for food and water consumption as well as for movement (Figure 18) with more feeding, 
drinking, and movement occurring during the dark cycle than during the light cycle. 

[0138] Comparison of each obese mutant group (2C, Ob, Ay) with their wild type 

controls revealed a significant effect of genotype for food consumption with all of the obese 
mutants consuming more food than their wild type controls (Figure 3). In addition, 
25 comparison of 2C mutants and WT mice also revealed a cycle by genotype interaction for 
food consumption, and surprisingly this interaction appears to be the result of the 2C 
mutants only consuming more food than their wild type controls during the light cycle 
(Figure 3). This is in contrast to the Ob and Ay mutants where the interactions of cycle and 
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genotype are not significant, indicating that the increased intake by the mutants relative to 
their wild type controls is similar in both the light and dark cycles. 

[0139] In contrast to food consumption, there was no significant effect of genotype 

for water consumption for any of the mutant groups compared with their wild type controls. 
5 However, 2C and Ob mutants exhibited significant cycle by genotype interactions. The 2C 
mutants exhibited a smaller decrease in water consumption from the dark cycle to the light 
cycle relative to their wild type controls. In contrast, the Ob mutants exhibited a larger 
decrease in water consumption from the dark cycle to the light cycle relative to their wild 
type controls. 

10 [0140] Examination of the quantity of movement revealed a significant effect of 

genotype for each mutant group compared with their wild type controls. The 2C mutants 
were hyperactive compared to their WT controls, and also exhibited a significant cycle by 
genotype interaction, indicating that the decrease in activity going from the dark cycle to the 
light cycle was greater for 2C mutants. Ob and Ay mutants were hypoactive compared to 

15 their respective WT controls, and both exhibited significant cycle by genotype interactions, 
revealing that the decrease in activity going from the dark to light cycle was less for the Ob 
and Ay mutants. 

Identification of Bouts and Clusters Using Feeding Behavior 
[0141] As seen above, there is clearly a diurnal pattern of ingestive behavior and 

20 physical activity that is altered in distinct ways by these different mutations. However, 

investigation of behavior on this time scale does not provide insights into the changes in the 
underlying behavioral units governing these alterations in diurnal intake. The examination 
of feeding behavior at a higher temporal resolution is necessary in order to identify the 
presence these behavioral units, bouts and clusters of feeding, and to determine their 

25 temporal pattern. In order to visualize such patterns, an event record of photobeam breaks, 
reflecting feeding, is displayed for a single 24-hour period for an animal of each genotype 
(Figure 19c). It is readily evident that the photobeam breaks are clustered together in time, 
and that the properties of these clusters vary both with time of day and with genotype. 

[0142] In order to quantify the presence of bouts and clusters of feeding, the gaps 

30 between photobeam breaks for individual animals were fit with log-normal distributions for 
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the light and dark cycles (Figure 20). It is clear from inspection of these frequency 
histograms that the gaps are best fit by more than two log-normal distributions in all cases. 
This indicates that there is a temporal structure underlying the generation of the gaps 
between photobeam breaks, and that the duration of these gaps is neither random nor 
5 regular. However, the pattern appears more complex than would be predicted if only bouts 
and clusters were present, since, if this were the case, only three log-normal distributions 
would be found: short intervals within bouts, longer intervals between bouts, and very long 
intervals separating clusters of bouts. Despite the presence of more than three distributions, 
the intervals have been classified into only three types. The reasoning behind this is that the 

10 first several distributions typically exhibit a very short duration, and the animal presumably 
remains in the feeder engaged in feeding behaviors rather than leaving the feeder and 
engaging other behaviors. This can be verified by determining the probability of the animal 
being in a position during these very short pauses, which have been labeled as within-bout 
intervals (WBIs) (Figure 21). This reveals that the probability of being at the feeder during 

15 the WBIs is much higher than the probability of being at any other location. It may be that 
the presence of several very short durations less than the bout criteria (DC (mean±sd): 
WT=0.4 min ± 0.2 min, =0.4 min ± 0.2 min, Ob=0A min ± 0.3 min, Ay =0.5 min ± 0.2 min; 
LC (mean±sd): WT=0.5 min ± 0.5 min, 2C=0.7 min ± 0.5 min, Ob=03 min ± 0.2 min, 
Ay=0.6 min ± 0.6 min) reflects different ways in which the animal handles the food while 

20 removing it from the feeder resulting in several different types of WBIs. 

[0143] Two additional distributions were routinely observed, inter-bout intervals 

(IBIs) and inter-cluster intervals (ICIs). IB Is are defined by being equal to or longer than 
the bout criteria but shorter than the cluster criteria, and ICIs are defined by being equal to 
or longer than the cluster criteria (DC: WT=19.5 min ± 7.9 min, 2C=19.4 min ± 6.3 min, 

25 Ofc=35.5 min ± 17.6 min, Ay=34.9 min ± 19.5 min; LC: WT=36.2 min ± 12.7 min, 2C=20.1 
min ±11.4 min, 06=36.1 min ± 21.6 min, Ay-A\.l min ± 14.9 min). If the IBIs and ICIs 
are truly distinct, then some feature of the animal's behavior, such as its location, should be 
different during the IBI as compared with the ICI. This is in fact appears to be the case as 
seen by examining the animal's probability of being in a position during the IBIs and ICIs. 

30 During the IBIs, the probability of the animal being around the feeder, around the licking 
spout, as well as elsewhere around the cage is roughly equal, and this distribution is 
markedly different than that observed for the WBIs. This is consistent with the idea that 
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during the IBIs the animal terminates a bout of feeding and engages in other behaviors. 
During the ICIs, the animal's probability of being at one location in the cage is much higher 
than being at any other location. This distribution appears distinct from that of the IBIs 
distribution indicating that during the ICIs and the IBIs the animal engages in distinct 
behaviors. In fact, daily recording of the animal's nest location reveals that the position of 
the nest is corresponds with the location in the cage that the animal has the highest 
probability of being at during the ICIs. This indicates that the bouts are clustered in time 
and are separated by longer intervals during which the animal is inactive at its nest. 

Comparison of Cluster Properties 
[0144] In comparison with their respective wild type controls, the obese mutant 

mice all exhibit increased food consumption, thus it is expected that the properties of the 
clusters during which the mice are consuming food should change in such a way as to 
account for this increased consumption. Thus during the light cycle, the 2C mutants exhibit 
an increase in the number of feeding clusters accompanied by a corresponding decrease in 
the duration of the ICIs, but with no significant difference in the cluster size compared with 
their wild type controls (Figure 22). In contrast, the Ob and the Ay mutants both show an 
increase in cluster size during the light cycle relative to their wild type controls. In addition, 
the Ob mutants also exhibit a decrease in the number of clusters accompanied by a 
corresponding increase in the duration of the ICIs. The Ay mutants show a similar trend of 
decreased cluster number and increased ICI duration relative to their controls but this does 
not reach significance. During the dark cycle, the Ob mutants again exhibit an increase in 
cluster size and a decrease in cluster number accompanied a increase in the duration of the 
ICIs. The Ay mutants show similar trends in the differences of their cluster properties 
relative to their WT controls however these do not reach significance. In contrast, during 
the dark cycle, the 2C mutants relative to their controls exhibit a somewhat distinct pattern 
of differences in their cluster properties in comparison with the light cycle differences. 
There is a significant decrease in the number of clusters but no significant change in cluster 
size or ICI duration. Examination of the circadian timing of certain cluster properties 
revealed that these properties are not constant across either the dark-light cycle nor within 
the dark cycle (Figure 23). WT mice as well as 2C and Ay mutant mice exhibit increased 
cluster sizes at the onset and offset of the dark cycle Ob mutant mice however have less 
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variation in cluster size within the dark cycle, seldom showing increased cluster sizes at the 
onset and offset of the dark cycle. All genotypes had relatively constant cluster sizes across 
the light cycle. Variations in cluster sizes within the dark cycle but not the light cycle may 
help explain why significant differences in cluster properties are more easily detected in the 
5 light cycle than the dark cycle. 

Quantification of Circadian and Ultradian Rhythms of Feeding Behavior 
[0145] The variation of cluster size within the dark cycle suggests that ultradian 

rhythms, periodicities shorter than 24 hours, are present in addition to circadian rhythms. 
Furthermore, it is apparent that differences in these circadian and ultradian rhythms are 
10 present across genotypes, making it important to quantify these periodicities allowing for 
comparisons between the genotypes. 

[0146] To quantify these periodicities, periodograms were generated from 10 days 

of data, and revealed significant peaks in amplitude not only at a period of 24 hours 
consistent with circadian rhythms but also at shorter durations consistent with ultradian 

15 rhythms in feeding (Figure 24). Notably the WT mice frequently exhibit ultradian peaks at 
12, 8 and 5 hours in addition to a prominent 24 hour peak while the Ob mutants have a less 
prominent 24 hour peak and exhibit fewer and shorter ultradian peaks. The significant 
periods were used to generate a waveform revealing that the onset and offset peaks in the 
WT appear to result from the combination of the 12, 8 and 5 hour peaks. Since these peaks 

20 are usually absent in the Ob mutants, the Ob mutants do not exhibit a strong onset-offset 
pattern of feeding and instead exhibt a more regular pattern of feeding at some interval 
dependent on the specific individual mouse. 

[0147] The periodograms were also used to cluster mice based on these patterns of 

feeding (Figure 25). Three qualitative groups were seen, one with only WT mice and 2C 

25 mutants, one with mostly WT mice and Ay mutants, and one with only Ob mutants. The 
waveforms generated from the averaged periodograms showed three distinct patterns of 
feeding. The WT and 2C mice showed higher onset than offset peaks in feeding. The WT 
and Ay mice showed higher offset than onset peaks in feeding. In contrast, the Ob mice 
showed a more periodic pattern of feeding with no distinct onset and offset peaks in 

30 feeding. These averaged waveforms show a pattern or feeding that is consistent with the 
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patterns of feeding that have been observed by examining the variations in cluster size with 
circadian time. 

Discussion 

[0148] The rationale for home cage behavioral monitoring stems from the 

5 hypothesis that behavior will be organized in time and space because such organization 
should improve the fitness of the animal. As a result, behavioral patterns, as with other 
phenotypes, will be genetically determined, and the description and quantification of these 
patterns should therefore lead to insights into the underlying function of the central nervous 
system structures regulating the production of these patterned behaviors. Here, data 

10 collection methods and analytical techniques for the quantification of behavioral patterns 
were developed and then tested by comparing animals with profound alterations in an 
important organizer of behavior, the maintenance of energy balance. Distinct alterations in 
patterns of ingestive behavior and physical activity were observed in mice having different 
mutations that result in obesity. These alterations can be understood in the context of 

15 genetic disruptions in a system attempting to maintain energy balance. 

[0149] A relatively simple model of ingestive behavior has been developed. It was 

observed that all mice organize episodes of feeding, and likely of drinking and movement, 
into bouts. The bouts themselves are clustered together in time, and these clusters are 
separated by prolonged periods of inactivity at the nest. In addition, the cluster size varies 
20 not only with a circadian but also with ultradian rhythms. There is thus both a temporal and 
related spatial structure to the organization of the behavior of a mouse in its home cage. 

[0150] The different mutant animals studied here exhibit distinct alterations in the 

underlying patterns that govern ingestive behavior. Thus, the 2C mutants were observed to 
have an increased number of clusters with a corresponding decrease in their ICI durations 

25 during the light cycle where they appear to be hyperphagic relative to their WT controls. 

However, the quantity of food that they consume during each cluster remains similar to that 
of their WT controls. This suggests that, during the light cycle, these mutant mice may a 
have disruption in systems regulating the initiation of clusters that results in clusters 
occurring more frequently with a shorter intervening interval. This might be interpreted as 

30 the result of an increased drive to consume food, but could also reflect disruptions of the 
sleep-wake cycle leading to increased arousal. This is interesting in light of the likely 
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connection between the regulation of energy balance and arousal suggested by the discovery 
of the orexin neuropeptides, which have been discovered to regulate both food consumption 
and patterns of sleep and arousal (Mieda and Yanagisawa (2002) Current Opinion in 
Neurobiology 12: 339-345). It should be noted however that while the increased frequency 
5 of clusters points to an alteration in initiation of behavior, especially an increase in the 

initiation offeeding, there is no compensating decrease in the amount of food consumed in 
each cluster which would adjust the overall intake during the light cycle to a level consistent 
with the WT controls. This might be because the 2C mutants also exhibit increased activity 
during the light cycle and thus require increased food intake as a result of their increased 

10 energy expenditure. However, the 2C mutants are also hyperactive relative to their WT 
controls during the dark cycle but do not exhibit increased intake during the dark cycle. 
This suggests that the effects of this mutation on food intake are dependent on the time 
point during the light-dark cycle and argues against a tight correlation between food intake 
and activity that would account for the increased energy intake observed in these mice. It 

15 may be that there are distinct metabolic processes that are selectively altered during the light 
cycle in these mice that lead to the eventual development of a mild obesity later in the life 
of these animals. Consistent with this, it is known that there are numerous alterations in the 
transcriptome of an animal across the light-dark cycle in organs, such as the liver, that play 
a major role in metabolic regulation (Akhtar et al. (2002) Current Biology 12: 540-550). In 

20 addition, it is also known that the 2C mutants display an age-related decrease in the cost of 
locomotion that is more marked WT mice and likely contributes to the development of 
obesity in the 2C mutants (Nonogaki et al. 2003). It would therefore likely be highly 
informative to examine how the patterns of ingestion and activity change in the 2C mutants 
with age. 

25 [0151] In contrast to the 2C mutants, the OB mutants exhibit a very distinct pattern 

of intake and activity relative to their WT controls. Both during the light and dark cycle, the 
OB mutants exhibit an increase in the quantity of food consumed during a cluster 
accompanied by a decrease in the number of clusters and a corresponding increase in the 
ICIs. It may therefore be the case that the OB mutants have an impairment of the system 

30 that regulates the termination of food intake. However, even though the frequency of their 
intake is decreased which would serve to decrease total intake, it is not sufficient to 
compensate for the increase in consumption during a cluster, indicating that the mechanisms 
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regulating the initiation of intake are also not functioning normally in this animal. Given 
the role of leptin in acting as a signal of adipose depot status, the concurrent alterations in 
the initiation and termination of intake may represent an appropriate CNS response 
attempting to increase energy intake. This is consistent with the profound hypoactivity 
5 observed in these animals as well as the numerous other metabolic and neuroendocrine 
alterations that they exhibit that seem to indicate a response to starvation (Ahima et al 
(1996) Nature 382: 250-252). Such a response to conserve energy however stands in 
distinct contrast to the patterns exhibited by the 2C mutants, who also develop obesity 
although later and of a much less pronounced degree. In addition, the OB mutants have 

10 another distinct alteration in their pattern of behavior that is not observed in either the 2C or 
the AY mutants, which is an alteration in both their pattern of circadian and ultradian 
rhythms. The OB mice have a less pronounced circadian rhythm, and in addition, do not 
exhibit the ultradian rhythms that give rise to the marked onset and offset peaks of feeding 
at the beginning and end of the dark cycle observed in the WT, 2C, and AY mice. Instead, 

15 the OB mutants tend to have a single shorter ultradian rhythm that give rise to a more 
regular periodic feeding pattern. One possible explanation for this pattern is that factors 
controlling energy balance dominate the regulation of the temporal patterns of activity and 
ingestion in the OB mutants overriding signals underlying the ultradian and circadian 
rhythms. The dampened circadian rhythm may reflect a decrease in responding to signals 

20 from the animal's internal circadian clock such as that located in the suprachiasmatic 

nucleus. The source of the signals generating the ultradian rhythms is not known but may 
also be overridden. Alternatively, it may be that the OB mutation is directly involved in 
altering the signals producing these periodicities. 

[0152] Finally, the AY mutants also exhibit distinct changes in their patterns of 

25 behaviors relative to their WT controls. These appear distinct from the patterns observed in 
the 2C and OB mice, although the AY mutants have some features of their behavioral 
pattern that are more similar to the OB than the 2C mutants. During the light cycle, the AY 
mutants have an increase in the amount of food consumed during a cluster that may account 
for their increased intake. However, no significant changes in cluster properties were 
30 observed during the dark cycle in the AY mutants although they appear to be as 

hyperphagic during the dark and light cycle compared to their WT controls. This lack of a 
detectable change in pattern during the dark cycle may be the result of the marked ultradian 
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rhythms (presence of onset and offest feeding peaks) present in both the WT and AY mice 
which makes comparison of cluster size across the dark cycle less meaningful. The 
presence of these onset-offset peaks in feeding during the dark cycle also seems to 
distinguish the AY pattern of ingestion from the OB pattern of ingestion. 

5 [0153] In fact, differences in the periodograms between mice that reflect these 

ultradian and circadian rhythms can be used to qualitatively group animals based on their 
periodicities of ingestive behavior. This reveals that the WT mice seem to be split into two 
groups: one with a more prominent onset feeding peak and another with a more prominent 
offset feeding peak occurring during the dark cycle. The 2C mutants appear to fall into this 
10 first group while the AY mutants fall into the second group, and for the most part the OB 
mutants fall into a third group which lacks clear onset and offset feeding peaks during the 
dark cycle. It should also be noted that the underlying cause of these onset and offset peaks 
appears to be the occurrence of long duration clusters that occur at the beginning and end of 
the dark cycle. 

15 [0154] In summary, a descriptive quantitative model of the home cage behavior of 

mice has begun to be developed. The properties of this model vary and exhibit distinct 
patterns as a result of different mutations producing obesity suggesting these mutations 
affect different systems regulating energy balance. In addition, these results indicate that it 
is possible to detect changes in the behavioral patterns of a freely behaving animal as it 

20 organizes its behavior in response to internal and external cues. The goal of such a 

description of behavioral patterns is to produce a model that reflects the structure of the 
data, which in this case is the continuous stream of behavior exhibited by an animal. 

[0155] In turn, the model hopefully reflects the structure of the generating process, 

here, the central nervous system. A compact model of behavioral data, such as the data 
25 collected and examined in this study, should allow us to learn about the functioning of the 
central nervous system and the underlying genetic contribution to its structure as well as to 
make predictions from the model that will be consistent with new samples. 

[0156] It is understood that the examples and embodiments described herein are for 

illustrative purposes only and that various modifications or changes in light thereof will be 
30 suggested to persons skilled in the art and are to be included within the spirit and purview of 
this application and scope of the appended claims. All publications, patents, and patent 

-46- 



305T-300600US 

applications cited herein are hereby incorporated by reference in their entirety for all 
purposes. 
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